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The respiratory chain or oxidative phosphorylation system (OxPhos) generates most of the chemical energy 
(ATP) used by our cells. The cytochrome c oxidase (CcO) is one of three protein complexes of OxPhos 
building up a proton gradient across the inner mitochondrial membrane, which is ultimately used by the 
ATP synthase to produce ATP. We present molecular dynamic simulations of CcO in a mimic of the 
mitochondrial membrane, and identify precise binding sites of cardiolipin (CL, signature phospholipid of 
mitochondria) on the protein surface. Two of these CL binding sites reveal pathways linking CLs to the 
entrance of the D and H proton channels across CcO. CLs being able to carry protons our results strongly 
support an involvement of CLs in the proton delivery machinery to CcO. The ubiquitous nature of CL 
interactions with the components of the OxPhos suggests that this delivery mechanism might extend to the 
other respiratory complexes. 



Mitochondria are the intra-cellular organelles that produce most of the chemical energy consumed by a cell 
and are therefore often refered to as the "power plant" of the cell. The energy is produced through the 
oxidative phosphorylation (OxPhos) system. The OxPhos system is embedded in the mitochondrial 
inner-membrane and consists of four large membrane protein assemblies, the so-called "respiratory chain 
complexes" (i.e.: NADH quinone oxidoreductase, complex I; cytochrome bc h complex III; cytochrome c oxidase, 
complex IV; ATP synthase, complex V) and by some small electron carriers (quinones and cytochrome c). 
Complexes I, III and IV operate a series of electron transfers whose tortuous paths through the respiratory chain 
system is coupled with the translocation of protons from the inside (matrix) to the mitochondrial intermembrane 
space (IMS), leading to an electro -chemical gradient ultimately utilized by the ATP synthase to transform ADP 
into ATP. 

Cardiolipin (CL), the signature phospholipid of mitochondrial membranes, has crucial implications in mito- 
chondrial processes, and on the OxPhos in particular 1,2 . CLs are formed by two phosphatidyl-glycerol molecules 
linked by an additional glycerol moiety and possess an anionic charge ( — 2 e). They are found in several parts of 
the mitochondria and constitute about 10% of the phospholipids of, for instance, the bovine heart mitochondrial 
inner membrane 3 ' 4 . CL deficiency correlates with numerous diseases 5 including the Barth syndrome 6 and heart 
failure 7 . Notably, the formation, stability, and function of individual respiratory chain complexes and of respir- 
atory chain supercomplexes strongly relies on the presence of CLs in the membrane 8 " 12 . CLs bound to individual 
complexes have been reported by several structural studies on various organisms 13 " 15 , leading to hypotheses on 
their possible involvement in a proton uptake pathway 16 and/or in assuring the structural integrity of individual 
complexes 16 " 19 and supercomplexes 15 . 

The focus of the present work is on bovine cytochrome c oxydase (CIV). It is a transmembrane protein complex 
that consists of 13 subunits labeled A to M as depicted in Fig. 1A (The subunit's nomenclature following a roman 
number code as defined by Kadenbach 20 is given in Supplementary Information). Primarily based on purification 
and delipidation by chromatography and binding affinity assays 131718,21 ' 22 two classes of CL binding sites have 
been defined: two sites with high-affinity and one to two additional sites with a lower affinity. The latter have been 
associated with the structural integrity of CIV and of its dimeric form because these CLs stabilize the subunits G 
and H (Via and VIb according to Kadenbach nomenclature 20 ), which are mandatory for the formation of the 
dimer 1718 . The two CLs binding CIV with high affinity have been associated with the regulation of the electron 
activity of the enzyme 21 ' 22 . 
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Figure 1 | Structural characteristics of bovine heart cytochrome c oxidase (CIV). (A) Structure of CIV with its thirteen subunits (A-M) shown in a 
cartoon representation with a chain-based color-code for one monomer. The subunit's nomenclature following a roman number code used in others 
studies 20 is given in Supplementary Information. We kept on the letter code to avoid confusion with the site names. In the bottom view of the dimer 
structure as found in the X-ray structure CLs are depicted 24 . There are three CL binding sites per monomer (CLla, CLlb and CL2), with one facing the 
matrix side (CL2) and two facing the intermembrane space (IMS; CLla and CL1 b). CL3 indicates the location of an additional binding site suggested by 
photolabeling experiments 19 . (B) Simulation box for CIV system, with the protein shown in green, POPC molecules in gray/white, CLs in red/orange and 
the aqueous phase in blue. (C) Structure of a cardiolipin molecule in an atomistic (AA) and a coarse-grain (CG) representation. 



It is only with the most recent crystal structures of CIV 23 ' 24 that CLs 
were found bound to the protein. In the dimeric form of CIV four 
CLs co -crystallize with the complex, defining three potential CL 
binding sites per complex (CLla, b and CL2 in Fig. 1A) after con- 
sideration of the dimer's symmetry. The intuitive assignment of the 
CL binding sites based on their location on the complex was in 
overall agreement with the in-depth photolabeling study of Sedlak 
et al. 19 . The low-affinity CL binding sites are located at the interface of 
the dimer and interact with subunits G (CLla/b in Fig. 1), while a 
high affinity binding site was found on a more membrane exposed 
protein surface interacting with subunit J (CL2 in Fig. 1). There is 
however a CL binding site identified by Sedlak et al. 19 that is missing 
in the crystal structure (CL3 in Fig. 1). 

Here we explore the binding of CL to complex IV in a monomeric 
form by means of coarse-grain molecular dynamics (CGMD) simu- 
lations. MD simulations have been successfully used to explore the 
lipid/protein interplay 25 ' 26 and more specifically lipid binding to a 
variety of membrane proteins 25 ' 27 " 32 . The use of a CG model presents 
the advantage to investigate the lipid binding to proteins on much 
longer time scales than conventional atomistic models and thus 
allows exploring the processes involved from a more realistic 
dynamic perspective. We recently applied a similar approach to 
the study of CL binding sites on the cytochrome bci 5 in which case 
we have characterized some key dynamic features of known CL 
binding sites and revealed a set of new membrane-exposed binding 
sites potentially involved in the formation of supercomplexes. 

The set of CGMD simulations presented here clearly identified the 
precise positions of the two CLs with high-affinity binding sites on 
CIV 19 . The lack of the low-affinity binding sites expected at the 
interface of the dimer are not observed in the simulations, which 



confirms that they are strongly associated with the dimer form of 
CIV and might only exist with the dimer. Five additional CL binding 
sites with low- affinity are found and may be easily rationalized in 
light of the position of other co- crystallized lipids 24 and their com- 
mon features. Most remarkably two of the binding sites are found 
located at the matrix entrance of known proton uptake pathways (D 
and H). We show how the structure of the protein reveals extensions 
of these pathways that directly link the CLs in these binding sites to 
the entrance of the pathways. In the context of the ability of CL to 
trap protons our data strongly support that CL maximizes CIV elec- 
tron transport activity by providing protons to the uptake pathways. 

Results 

Identification of seven CL binding sites on bovine heart cyto- 
chrome c oxidase. We characterized the CL density around the 
protein from a 100 jas CGMD simulation of this complex embed- 
ded in a POPC/CL (20/1 molar ratio) membrane. The average 
densities of CLs (Fig. 2 A) demonstrate the existence of several 
preferential sites of interaction of CL with the protein. We define 
CL binding sites as the locations having a CL density more than five 
times the bulk value. Accordingly, seven binding sites are found on 
either leaflet of the membrane and are labeled I to VII (Fig. 2A). Sites 
I-V are on the matrix side of the membrane and sites VI and VII on 
the side of the IMS. Sites V and VII may be occupied by two CLs and 
are therefore subdivided into Va/Vb and Vlla/VIIb. A hierarchy of 
binding strength was observed, starting with many CL binding/ 
unbinding events for some sites (e.g. 380 exchanges for site Vllb) 
but only few exchanges for others (e.g. few unbinding events for site 
I) during the same amount of simulation time. For all sites but site I, 
binding events by at least 4 CLs were observed (see Supplementary 
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Figure 2 | Cardiolipin (CL) binding sites on cytochrome c oxidase (CIV). Binding sites are extracted from a 100 [is of CGMD simulation of the 
complex embedded in a CL/POPC membrane bilayer. (A) From left to right: matrix, membrane (two orientations) and inter membrane space (IMS) 
views of CIV with the CL densities shown in yellow volume maps at an isovalue corresponding to at least 5 times the bulk density. The protein is shown as 
shaded grey cylinders with the CL densities projected onto them. (B-H) Detailed description of the CL binding sites I to VII, respectively. The residues are 
numbered as follows: "chain:residue sub Slte ". For each site, the subunits involved in the interactions with the CLs are depicted as colored cartoon as in 
Fig. 1. The rest of the protein is shown in a transparent gray cartoon. 



Table 1 | Occupation (H) and residence time (9, (is) of cardiolipin binding sites on cytochrome c oxydase (CIV). The values are averaged 
over a 1 00 jis CGMD simulation. The accuracy of the occupation levels is ±0.02 at most and of the residence times within ±0.1 jis. See 
the Extended Methods given as Supplementary Information for details. #CLs is the number of different CLs getting in contact with a site 
and #events the number of binding/unbinding events. The lipids described by Shinzawa-ltoh et ol. 24 (see Fig. 4) at the locations 
corresponding to the CL binding sites found in the simulations are indicated in the bottom row. The CL binding site predicted by 
Sedlak efa/. 19 at the site II is also indicated 
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CL 


TGL2 


PE3 


PE1 


TGL3 


TGL3 




TGL1 


TGL1 



The strong binding of CLs in sites I and II led to poor binding/ unbinding events statistics and prevented us from determining accurate residence time. Instead we indicate a rough estimate obtained from the 
binding behavior. See Figure SI in Supplementary Information. 
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Movie), allowing us to extract meaningful statistical averaged 
occupations and CL residence times listed in Table 1. The detailed 
structural characterization of the binding sites is presented in 
Figure 2B-H. 

Site I quickly bound a single CL that remained bound during the 
entire simulation although reorientation of the CL led to multiple 
short unbinding events (see Supplementary Figure SI). These events 
complicated the computation of the lifetime estimated at more than 
50 |is (Table 1). In this site the CL mainly interacts with the subunit C 
and with a couple of residues of subunit J (Fig. 2B). Sites II and III are 
adjacent to site I and located at the junction of subunits A and L and 
subunits C and G, respectively (Fig. 2C-D). They both bind one CL at 
a time and have slightly lower occupancies (60%) than site I. 
Although it is difficult to precisely determine the time scale of CLs 
binding to site II (see Supplementary Fig. SI; estimate —60 jas) it is 
significantly longer than in site III in which CLs exchange on a time 
scale of ~ 1 |is. Site IV is located on the face of CIV that corresponds 
to the dimer interface in the crystal structure. This site, involving 
subunits A and C, is occupied only 36% of the time and CLs bind on a 
sub -microsecond timescale. Site V is located at the other extremity of 
CIV as compared to site I. CLs in site Va and Vb interact with 
subunits B, E and I, and A and B, respectively. They have similar 
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Figure 3 | Potential of mean force for binding of various lipids to sites I, II 

and III. (A) Comparison of CL's binding strength to site I (cyan), II (blue) 
and III (green). (B) Comparison of binding strength of TGL (yellow), 
POPC (orange) and CL to site II. Two CLs were tested; double charged 
(— 2 e) and neutral, blue and green curves, respectively. In both panel the 
relative free energy of the system is expressed as a function of the distance, 
dcoM> between the center of masse of the lipid headgroup and of the 
binding site as defined in Fig. 2. The error on the measure (estimated using 
the Bayesian bootstrapping method) is shown by the shaded area behind 
the curves. 



occupancies, —70%, and similar residence times, —0.5 jas. Site VI is 
located on the IMS side of CIV. It involves the C-termini of subunits 
L and M, has a low CL occupancy, 50%, and a small CL residence 
time, —0.4 |is. Site VII is located on the same protein surface as site 
V, but on the opposite leaflet (IMS). CLs in site Vila and Vllb 
primarily interact with subunits A, B, D and I on the IMS side of 
the protein. Site Vila is significantly more occupied (88%) than site 
Vllb (50%) and has a relatively long residence time of —10 jas. The 
difference between sites Vila and Vllb might reflect that site Vila is 
buried deeper in a cavity formed by helices of subunits B and I, 
whereas site Vllb is more exposed to the membrane. 

We further characterized the binding strength of CL to CIV by 
determining the potentials of mean force (PMFs) of CL's binding to 
a selection of the sites (Fig. 3). The data clearly indicates that CL's 
binding to site I and II is stronger than to the other sites. Site III was 
taken as a reference site because it is relatively well occupied, accessible 
to the membrane and explored by many CLs so that the statistics 
obtained is reliable (Table 1). The data also indicate that CL binds 
similarly to site II than to site I. The comparison of the PMFs of a CL 
binding to site II with or without the — 2 e charge of its headgroup 
indicates that although the charge of the headgroup bears the most 
significant contribution to the binding, the tails also contribute to CL's 
binding strength. The contribution from the CL's bulky tail is by itself 
larger than the binding strength of a POPC molecule. 

Residue content of the binding sites. Interesting features emerged 
from the analysis of the residues in contact with CLs during the 100 
|is simulation (Fig. 4). First, all the residues that made at least one 
contact with a CL molecule virtually span the entire transmem- 
brane region of the protein, indicating that the CLs explored the 
complete transmembrane protein surface. As expected, the distri- 
bution of contacts as a function of the residue type (grey bars in 
Fig. 4) is in general agreement with the amino acid distribution in 
integral membrane proteins 33 keeping in mind that a residue would 
need to be exposed to the membrane to be in contact with a CL. 
Furthermore, focusing on the composition of the CL binding sites 
(black bars in Fig. 4), the positively charged amino acids (lysine and 
arginine) were found to account for —25% of the CL ligands. The 
large contribution from positively charged residues might be 
expected since CL carries a -2e charge and is illustrated by the 
strong correlation between the locations of the CL's binding sites 
with the positive regions of electrostatic potential on the protein's 



60 



50 



g 40 

CD 

o 

§ 30 



o 
o 

O 20 



10 




llln 



.I.I 



M. 



limn.. 



Lys Leu His Tyr Trp Asp lie Asn Ser Cys 
Arg Phe Met Gin Ala Gly Thr Pro Val Glu 

Figure 4 | Residue content of the CL binding sites of cytochrome c 
oxidase (CIV). The gray sticks indicate the percentage of each residue type 
at least once in contact with a CL; e.g. present in the section of the CIV 
accessible to CLs. The black sticks give for each residue type the percentage 
of its participation to the CL binding sites. 
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surface (Fig. S2). Lysine is slightly favored over arginine (14 vs 11%), 
which contrasts with the results obtained previously for Oil 15 . In 
that case arginine was significantly more represented in the binding 
sites than lysine, which then represented only —5% of the CL's 
ligands against 18% for arginine. Phenylalanine and leucine were 
found to be a relatively important contribution to the CL binding 
sites (> 10%) in both CIII and CIV. 

Discussion 

A thorough analysis of the lipids bound to CIV was presented by 
Shinzawa-Itoh et al. 24 combining a crystal structure of the complex 
with mass spectroscopy. They discussed up to seven species of lipids 
per monomer of CIV among which three CLs, three phosphatidy- 
lethanolamines (PEs), four phosphatidylglycerols (PGs) and three 
triglycerides (TGLs). These lipid types are found in the native mito- 
chondrial inner membrane 3 (TGL only recently 24 ) and their presence 
in the crystal is therefore not expected to be an artifact of the crys- 
tallization process. Their positions are depicted in Figure 5 together 
with the CL binding sites found in the simulations. The comparison 
of the two sets of molecules allows to unambiguously identifying the 
two CL binding sites with high-affinity and involved in CIV's proton 
transport activity, and suggesting the most likely candidates for the 
low-affinity binding sites that have been linked to CIV's structural 
integrity. Our interpretation is in complete agreement with earlier 
predictions of Robinson and co-workers 13,19 ' 22 . 

Among the three CL binding sites per monomer found in the 
dimer structure of CIV 23 ' 24 (CLla, b and CL2 in Fig. 5), only CL2 is 
observed in the simulations. It is located at the site I. This site is filled 
early in our simulation and remains occupied during the entire 
simulation (Table 1 and Supplementary Fig. SI), which gives it one 
of the highest binding affinities observed in the simulations. This 
result is in agreement with its presence in the crystal structure and 



Matrix view IMS view 




Figure 5 | Location of the co-crystallized phospholipid positions and 
densities computed with our CGMD simulation. The experimental 
positions (upper part) were extracted from the PBD entry 2dyr 24 . The 
densities (lower part) are extracted from the 100 (as CGMD simulation. 
Views from the matrix (left) and the IMS (right) sides of CIV are shown. 
CL: cardiolipin; PC: phosphatidylcholine; PE: phosphatidylethanolamine; 
PG: phosphatidylglycerol; TGL: triglyceride. 



the PMFs described above, and confirms that it corresponds to one of 
the high-affinity sites as predicted by Sedlak et al 19 . CLla/b, located 
at the dimer interface in the crystal structures 23 ' 24 , are not stable in 
our simulations of CIV as a monomer. These two CLs have been 
suggested to be the ones loosely bound to CIV and strongly con- 
nected to the existence of the dimeric form of CIV by their stabiliza- 
tion of subunits G and H 1819 . Our simulations indicate that the tight 
association of CIVs in the dimer is a determinant factor for the 
stability of CLla and CLlb. To test this hypothesis, additional simu- 
lation of the dimeric structure of CIV including the experimental 
CLla/CLlb were performed. Both sites were found stable on the |is 
time scale. 

The site II defined by the simulations (Fig. 2 and 5) confirms the 
location of the second high- affinity CL binding site predicted by 
Sedlak et al. 19 using photolabeling and missing in the crystal struc- 
tures (CL3 in Fig. 1). The relatively low occupancy and residence 
time (Table 1) might be surprising at first but reflect that it took 
—40 |is for a CL to enter the protein cavity forming the binding site 
(see Supplementary Fig. SI). Once the CL made the proper contacts 
the site is fully occupied and therefore corresponds to a high-affinity 
binding site. This result is corroborated by the PMF of CL binding to 
site II, which shows a similar bonding strength to site I (Fig. 3A). It is 
interesting to note that a TGL (TGL2 in Fig. 5) occupies the site II in 
the crystal structure 24 . From our understanding of the system it is 
extremely unlikely that both CL and TGL molecules occupy site II. 
The PMFs of the binding of a CL vs. a TGL (Fig. 3B) demonstrate not 
only that CL binds much stronger but also that TGL is no stable in 
site II. Unrestrained simulations confirm this observation: a TGL 
placed in site II left it quickly and went to mix with the bulk after 
exploring the surrounding of the protein surface around the site. The 
PMFs showed that an important part of the CL's binding strength is 
due to coulombic interactions, which might rationalize the weak 
binding of TGL; it is a neutral molecule. It is not clear at this point 
why a TGL molecule is present in the crystal although it was argued 
not to be an artifact of the purification 24 . One might speculate that the 
TGL molecule had the capability to replace the CL using the similar- 
ity of their bulky tails. This scenario implies however that the CL was 
first destabilized, which might occur by disruption of the interaction 
of its headgroup during the crystallization process. Note the PMFs 
showed a neutral CL still binds site II stronger than TGL suggesting 
that TGL's bulky tails are not sufficient to stabilize it. See the 
Supplementary Information for additional discussion on the pres- 
ence of TGL molecules in CL's binding sites. 

The remaining CL binding sites (III-VII) found in the simulations 
are relatively weak binders and unless they are buried within a large 
protein cavity (sites Va and Vila) they have a low occupancy. They all 
correspond to the location of binding of another lipid type in the 
crystal, which is in line with previous studies showing the conser- 
vation of binding sites by different lipid types 14 . We review the several 
cases observed in CIV and show that it is quite straightforward to 
rationalize the discrepancies from the similarities of CLs with PGs 
and TGLs (CL and PG share a negatively charged headgroup, while 
CL and TGL bulky tails), see Supplementary Information for more 
detailed discussion. 

The high degree of conservation of lipid binding sites on CIV from 
various organisms 14 suggests an important functional role of these 
bound lipids. CLs have been shown to affect CIV in two ways. First, 
there is extensive data demonstrating that the presence of two CLs 
with high-affinity binding to CIV, which we have shown are bound in 
site I and II, is mandatory for a fully functional CIV 22 . It is however 
unknown by which mechanism they operate. CIV enzymatic activity 
is often summarized by its electron transport activity (up take from 
cytochrome c) but CIV uses as many protons as electrons to trans- 
form a dioxygen molecule into two water molecules. CLs have the 
ability to trap protons 34 and thereby can facilitate proton transloca- 
tion along the membrane surface 35 . This particular feature of CL has 
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been suggested to be an important aspect of Oil mechanism for 
proton transfer 16,36,37 and might extend to CIV by providing proton 
sources to either of the D-, K- or H-pathways 38 " 41 . The proximity of 
CL3 to the D -pathway entrance has suggested that CL might be 
acting as a proton antenna to this pathway 19 although only limited 
structural insight is available since this site is not occupied in the 
crystal structures 23 ' 24 . 

Our data confirm that this ability of CL to carry protons is indeed 
relevant to site II but might also be to site Vb. Both sites are located on 
the matrix side of the IM where the protons are taken up. Site II 
corresponds to CL3 (Fig. 1, 5) and is very close to the D-pathway 
entrance (—1.1 nm; Fig. 6), while site Vb is close to the entrance of 
the H-pathway (—0.8 nm; Fig. 6). Up to now residues A:D91 
(A:D132 in Rodobacter sphaeroides) and A:D407 define the matrix 
side entrances to the D- and H-pathways, respectively. Only little is 
known on the proton path before these residues. A close inspection of 
the protein structure revealed in both cases the existence of a strong 
H-bond network starting at the CL binding sites II and Vb and 
leading to the entrances of the D- and H-pathways, respectively 
(Fig. 6). These networks form clear extensions of the D- and H- 
pathways towards the exterior of the protein on the matrix side. 
They directly connect the CL binding sites to A:D91 (A:D132 in 
Rodobacter sphaeroides) and A:D407, respectively, and thereby 
strongly support the idea that CLs provide a source of protons at 
the surface of the membrane facilitating CIV electron transport 
activity. 

It is also notable that up to four CLs occupy the protein cavity close 
to the H-pathway; two lipids in sites Va/b and Vlla/b on the matrix 
and IMS side of the membrane, respectively. This aggregation of CLs 
might be relevant to the 'dielectric channel' activity of CIV proposed 
by Rich 42 . 



The second way by which CLs affect CIV is by providing structural 
integrity. Two additional CLs, which we propose being CLla and 
CLlfr in the crystal structures (Fig. 1,5), serve the structural integrity 
of CIV by stabilizing subunits G and H, and thereby its dimeric 
form 17 . These two CLs have been shown to bind CIV with a low- 
affmity 19 and are not observed in the simulation of the monomeric 
form of CIV (Fig. 2) but are stable in the dimeric structure. This 
ability of CLs to stabilize CIV dimers can be extended to high-order 
oligomers of the respiratory chain complexes as it was reported in the 
context of the formation of supercomplexes. This behavior is most 
notable for the formation of supercomplexes CIII 2 -CIV and CIII 2 - 
CIV 2 912 . Our recent CGMD simulations of the CIII 15 suggested that 
CL binding sites on the surface of the complex might define the 
location of protein-protein contact with CIV. The CL binding sites 
found on CIV might share a similar function. 

In summary, the CGMD simulations used in this study to invest- 
igate the CL binding sites on the respiratory chain complex IV, 
cytochrome c oxidase, provide new and valuable insights on the 
way CLs participate to the function of proteins. The locations of 
the CL binding sites that most likely correspond to both the high 
and low- affinity CL binding to CIV are identified with high level of 
confidence. They agree with and reconcile all known experimental 
data. CL binding sites on the surface of the protein are found at the 
proximity of two of the three known proton-uptake pathways, the D- 
and H-pathways. Two clear interaction networks connecting the CL 
binding sites to the entrances of these two pathways are uncovered. 
To our knowledge they provide the first complete proton pathway 
between the membrane surface (CL binding sites) and the D- and H- 
pathways in cytochrome c oxidase. This data strongly supports that 
CLs take active part in the proton delivery mechanism. Given the 
wide impact of CL in the respiratory chain machinery the role 



D-pathway 



H-pathway 



A:W6 
A:R5 



A:D14 



A:M339 b 
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Figure 6 | H-bond networks leading from the cardiolipin (CL) binding sites to the entrance of D- and H-pathways, A:D91 and A:D407, respectively. 

The CL binding sites are shown in stick representation. The yellow surface maps depict the CL densities shown in Fig. 2. A:D91 and A:D407 are shown in 
large sticks. The residues and water molecules (red spheres) participating to the networks are shown in a ball- and- stick representation and numbered as in 
Fig. 2. The bottom row shows side and top views of CIV with the CL densities, the residues involved in the sites in blue, the residues involved in the 
transmembrane section of the proton pathways (red arrows) in orange spheres and the heme molecules in green. The large red spheres position the 
entrances (A:D91 andA:D407) of the pathways. 
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suggested by our data should extend to other components of the 
respiratory chain. 

Methods 

Molecular models. The Protein Data Bank (PDB) entries locc and 2occ 38 were used 
to built a complete atomistic model of CIV, which was simulated as a monomer, its 
functional form. In the experimental dimeric structure of CIV 24 four CLs co- 
crystallize with the protein (Fig. 1 A). In a monomeric configuration, the experimental 
CL positions are all exposed to the membrane bulk and were not included in the initial 
conformation. In the simulations, the complex lipid composition of the 
mitochondrial membrane was modeled by a two components mixture of CL:POPC 
with a molecular ratio of 1:20 (5% of CL, 10% of total phosphorus content). This value 
is close to the experimental molecular ratio observed for bovine heart mitochondria 
(15 to 20% of the phosphorus content 3 ). 

The systems were described with the Martini CG force field for biomolecules 
(version 2.0 43 ) and its extension to proteins (version 2.1 44 ) together with the EINeDyn 
approach 45 . The Martini force field defines chemical groups as CG beads (or super- 
atoms) parameterized to reproduce known thermodynamical observables associated 
to these groups. Molecules are constructed from these super-atoms using conven- 
tional bond, angle and dihedral potentials such that they reproduce the structure and 
flexibility of atomistic models 43 ' 44 . EINeDyn maintains the secondary and tertiary 
structures of CG proteins through an elastic network. In our model, each subunit is 
maintained independently, so that domain motions (reorientation of subunits) are 
possible. CG parameters for CLs and TGLs were taken from the works of Dahlberg 
et al. 46 and Vuorela et al. 47 . 

Simulation details. All CG simulations were performed using the GROMACS 
simulation package version 4.0 48 . Conventional simulation setups associated with the 
use of Martini and EINeDyn were used. The protein, membrane bilayer (POPC and 
CL) and aqueous phase (water and Na + ions) were coupled independently to an 
external temperature bath 49 at 300 K and the pressure was weakly coupled to an 
external bath 49 at 1 bar using a semi-isotropic pressure scheme. Further details of the 
models, simulation protocols and limitations, as well as the methods used for analysis, 
are published as Supplementary Information. 

The CG protein structure was inserted into a pre-equilibrated CL:POPC mem- 
brane patch. The system was energy minimized and simulated for 10 ns with position 
restraints on the backbone beads of the protein to relax the solvent and side-chains 
before starting production runs. The system shown in Fig. IB contains the protein 
(1780 residues; 4117 beads), a POPC bilayer (966 lipids; 12,558 beads) including CLs 
(48 CLs; 1296 beads) and the aqueous phase (51,549 water beads and 103 sodium 
ions). 

The potentials of mean force (PMF) of lipid molecules binding to CIV binding sites 
were computed using an umbrella sampling approach with the distance to the protein 
surface as reaction coordinate. Details are given as Supplementary Information. 
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